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ABSTRACT 

We analyze and model an M1.0 flare observed by SDO/AIA and Hinode/ElS 
to investigate how flare loops are heated and evolve subsequently. The flare is 
composed of two distinctive loop systems observed in EUV images. The UV 1600 
A emission at the feet of these loops exhibits a rapid rise, followed by enhanced 
emission in different EUV channels observed by AIA and EIS. Such behavior is 
indicative of impulsive energy deposit and the subsequent response in overlying 
coronal loops that evolve through different temperatures. Using the method we 
recently developed, we infer empirical heating functions from the rapid rise of 
the UV light curves for the two loop systems, respectively, treating them as two 
big loops of cross-sectional area 5" by 5", and compute th e plasma evolution 



in the loops using the EBTEL model (IKlimchuk et al.ll2008l ). We compute the 



synthetic EUV light curves, which, with the limitation of the model, reasonably 
agree with observed light curves obtained in multiple AIA channels and EIS 
lines: they show the same evolution trend and their magnitudes are comparable 
by within a factor of two. Furthermore, we also compare the computed mean 
enthalpy flow velocity with the Doppler shift measurements by EIS during the 
decay phase of the two loops. Our results suggest that the two different loops 
with different heating functions as inferred from their footpoint UV emission, 
combined with their different lengths as measured from imaging observations, 
give rise to different coronal plasma evolution patterns captured both in the 
model and observations. 



Subject headings: Sun: corona - Sun: flare - Sun: magnetic reconnection - Sun: 
UV radiation 
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Introduction 



It is widely accepted that magneti c reconnection occurs in solar flares with a large 
amount of energy released in the corona (IPriest &: Forbes! 120021 ). By reconnection, magnetic 
field lines change connectivity and form new loops. The released energy is transported down- 
ward along the loop by non-thermal electrons or thermal conduction to produce enhanced 
optical and ultraviolet (UV) emissions at the footpoints of the loop. In the meanwhile, 
impulsive energy deposition in the lower atmosphere evaporates hot plasma into the newly 
formed flare loops. The plasma then cools down, becoming visible sequentially as loops in 
soft X-ray and extreme-UV (EUV) emissions. 

Assuming that energy deposit takes place during the impulsive phase of the flare, some 
modelers adopt a loop heating rate proportional to the observed har d X-ray light curve and 



model heating and evolution of the flare as one single loop (e.g., iFisher &: Hawleyl 11990 
Raftery et al.l 120091 ). However, it has b een recogn i zed t h at most flares consist of multiple 



loops formed and heate d successively ( IHori et al. 



19981: lAschwanden fc Nightingale] 12005 



Mulu- Moore et al.ll201ll ). With this consideration, IWarren ( 20061 ) applied multi-thread ID 
hydrodynamic simulations to a flare, and was able to reproduce the soft X-ray radiation 
of the flare as observed by GOES and Yohkoh. To match observations, the author used 
50 threads (flare loops) , each heat e d at a different time for about 200 s by a prescribed 
heating rate. Similarly, Hock et al. ( 20121) have modeled a two-phase flare observed by the 



Atmos pheric Imaging Ass embly (AIA; iLemen et al.ll2012 ) and EUV Variability Experiment 
(EVE: I Woods et al.ll2012[ ) on board the Solar Dynamics Observatory (SDO). They observed 



different flare loop systems brightened in these two pha ses, and used the 0D model called 
"enthalpy-based thermal evolution of loops" (EBTEL; iKlimchuk et al.l 120081 ) to compute 
plasma evolution in these loops with arbitrary heating rates that are adjusted to match with 
the observed EUV radiation. In their study, 22 loops are heated throughout the flare to 
match the observed flare evolution. 

Key to all loop heating models is the heating function, which describes where, when, 
for how long, and by how much a flare loop is heated. To date, most multi-thread loop 
heating models prescribe an arbitrary number of loops, each given an arbitrary amount of 
heating at an arbitrary time, and the parameters are adjusted to best match the observed 
light curves of coronal emission. To relate fla re heating with magnetic energy release in a 
more realistic manner, lLongcope et al.l (120 lOf ) developed a model to compute the rate of 
energy released in reconnection using the time-dependent reconnection flux measured from 
observations of flare ribbon evolution. The amount of released energy is then used to compute 
flare plasma evoluti on and coronal r adiation and compare with observations by GOES and 
RHESSI. Recently, iQiu et al.l ( 120121 ) introduced an alternative method to observationally 
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infer heating rates by analyzing spatially resolved UV light curves emitted at the feet of 
over a thousand flare loops anchored to brightened UV footpoints. The loop heating occurs 
when the UV emission starts to rise, the duration of heating is twice the rise time of the UV 
light curve, and the magnitude of the heating rate is assumed to be directly proportional to 
the UV counts measured at the foot of the loop. With this method, many properties of the 
heating rate are observationally constrained. 

For its high efficiency in modeling a large number of loops, a number of studies have used 
the EBTEL model to compute plasma evolution in flare loops, as listed in Table [TJ These 
studies by different groups treat the model inputs and outputs in different ways. The input 
to the model includes the number of flare loops and their properties (such as the loop length) 
and their heating functions. The model-calculated temperature and density of flare loops are 
then used to compute X-ray and EUV emissions by flare loops that can be compared with 
observations. Furthermore, the EBTEL model includes a mean enthalpy flux between the 
corona and its base at the transition region, which describes upflows predominant during the 
heating phase and downflows during the cooling phase. This is consistent with the picture 



phase of the flare, as repon 


ted in observations ( 


Antonucci & Dennis! 


19831; lAntonucci et al. 


1985. 1990: Schmieder e 


t al. 


1987; 


Ding et a 


1996: 


Brosiuf 


J2003|: 


Rafterv et al. 


20091) and ID 


hydro dynamic models ( 


Doschek et al. 


1982; 


Fisher et al. 


1985a 


b). Therefore, the EBTEL 



computed flow can be also compared with Doppler shift measurements to provide another 
constraint to flare heating models. 



Table 1: Comparison among the studies referring to the EBTEL model and flare 

HEATING 



heating function 



loop number 



compare with observations 
(instruments) 



Rafterv et al. (2009 



inferred from hard X-ray light curve 



temperature, emission measure, velocity 
(RHESSI, GOES, TRACE, CDS) 



Hock et al. (2012 



free parameters 



EUV light curves 
(EVE) 



'iu et al. (2012) 



inferred from spatially resolved UV light curves 



X-ray and EUV light curves 
{GOES, AIA) 



this study 



inferred from spatially resolved UV light curves 



EUV light curves, velocity 
(AIA, EIS) 



Following iQiu et al.l (120121 ). we will observationally infer heating functions from UV 
light curves and compute plasma evolution in individual flare loops. It, however, takes 
a few steps ahead of previous studies by making use of the most stringent observational 
constraints, in both the input heating functions and the output radiation and flow signatures. 
The novel aspects in this study include (1) the heating functions and predicted coronal 
radiations are derived and compared for different loops that are identified from imaging 
EUV observations and exhibit distinctive evolution patterns; (2) both radiation and velocity 
signatures are compared with observations in multiple lines by multiple instruments, AIA 
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and the EUV Imaging Spectrometer (EIS; ICulhane et al.l 120071 ) on board Hinode. Such 
more comprehensive analysis and modeling allow us to bette r constrain the lo op heating 
model and diagnose flare heating more accurately. Whereas iQiu et al.l (120 121 ) have used 
EBTEL to model a few thousand loops with the same set of free parameters, in this study, 
we will focus on modeling and analyzing only two loops in order to refine determination 
of heating functions. The 0D EBTEL allows us to investigate more efficiently variations 
of these free parameters and gain insight into the underlying physics that governs these 
parameters. Among these parameters, one is related to the heating rate, and the other 
describes the energy loss through the lower atmosphere. These parameters might not be the 
same in different loop systems, and will be examined in the present study. It also should be 
noted that the 0D EBTEL code is benc h-marked by more sop histicated ID hydro dynamical 
simulations with reasonable agreement (IKlimchuk et al.l 120081 ); on the other hand, with the 
EBTEL model combined with our method, we will be able to determine the first-order 
heating rates to be used as a valuable reference for more sophisticated ID hydrodynamic 
models. In the following Section [2j we describe flare observations. The model computation 
and comparison with observations are presented in Section |3j We give the summary and 
discussions in Section [U 



2. Observations and Data Reduction 

Observations of a GOES M1.0 class flare on 2011 February 16 are summarized in Figure 
[TJ The top panel shows normalized light curves of the flare in GOES 1-8 A, UV 1600 A, 
and EUV 193 A bands. UV and EUV images in the middle panels were obtained by AIA 
with the pixel scale of 0.6" and time cadences of 24 s and 12 s, respectively. Seen from 
the figure, the flare exhibits a rapid rise in the UV 1600 A. The UV emission is from the 
upper- chromosphere or transition region, where flare loops are anchored. It is followed by 
GOES soft X-ray emission, and then the EUV emission at 193 A. This pattern indicates 
an impulsive heating of flare plasma and subsequent cooling toward lower temperatures. 
From AIA imaging observations, we can identify flare loops and their footpoints to study 
the heating and evolution of individual flare loops. 

The flare was also observed by EIS. The 2" slit of EIS scanned the flaring region from 
west to east using a 5" step mode with an interval of 5 minutes and 50 seconds. The time 
ranges of five scans during the flare are indicated by the vertical dashed lines in the top panel 
of Figure [TJ These scans cover a field of view of 180" x 160" as indicated by the black box in 
the middle-left image in Figure [TJ The EIS data is reduced using the standard EIS software 
data reduction package. The software corrects detector bias and dark current, as well as hot 
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pixels and cosmic ray hits, and yields absolute intensities in units of erg cm -2 s _1 sr _1 A -1 . 
In this study, we focus on analyzing the Fe XII 195.12 line (formed at ~1.25 MK) and the Fe 
XV 284.16 line (~2.0 MK) because of the high signal-to-noise ratio in these two lines. When 
measuring the Doppler velocity, we use the average line center over the EIS observing region 
as the reference wavelength. 

We co-align images from different instruments. The AIA image is co-aligned with the 
HMI magnetogram using the aia_prep subroutine. The co-alignment between AIA and EIS 
observations is made using the 193 A image observed by AIA and 195 A image observed by 
EIS. 



3. Heating and Evolution of Two Flare Loops 

High resolution AIA imaging observations allow us to identify flare loops and their 
footpoints, and study their evolution from the lower-atmosphere to corona. The morphology 
evolution suggests that the flare is composed of two loop systems that are formed and heated 
sequentially. In this study, we treat these as two big loops, and identify four footpoints where 
these two loops are rooted. These four footpoints, named hereafter as FP1, FP2, FP3, and 
FP4, are indicated in the middle- left panel in Figure 1. They each have an area of 5" x 5" 
(limited by the EIS spatial resolution of 5" x 1"). FP1 and FP2 are located in the positive 
magnetic field region, and FP3 and FP4 are in the negative field region. Clearly shown in 
the EUV images (the middle-right panel in Figure 1), FP2 and FP3 are connected by a 
relatively short loop (called LP23 hereafter), and FP1 and FP4 by a longer loop (LP14). In 
the bottom panels of Figure [T] we show the EIS intensity and Doppler velocity maps in the 
Fe XV line. The four footpoints (and the two flaring loops) are also indicated in these maps. 

To study heating and evolution of these two loops, we adopt our recently developed 
method to infer heating rates from the UV 1600 A emission at the four footpoints, and then 
use a 0D loop heating model to compute plasma evolution in the two loops. 

3.1. Evolution of Flare Plasmas with Observationally Inferred Heating 

Functions 



A 0D EBTEL model (IKlimchuk et al.ll2008l ) is used to compute the coronal response to 
a specified energy input. The EBTEL model solves time-dependent energy and momentum 
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equations to compute the mean temperature T and density n of coronal plasmas: 

dn 2 c 2 

dt 5kBT 



§ + cm 2 A(T) 



[Q(t)-(l + Cl )n 2 A(T)} 



(2) 



dp _ 2 
~dl ~ 3 

where p = 2 ksfiT is the pressure, fee is Boltzmann's constant, C2=0.87 is the typical ratio 
of mean temperature to apex temperature, L is the loop half-length (note that the model 
assumes a symmetric loop, so only the half-length is considered), Q(t) is the volumetric heat- 
ing rate, A(T) is the radiative loss function, and F c is the thermal conductive flux. Thermal 
conduction is computed with a classical Spitzer-Harm form, F c = — | k (T/C2) 7 / 2 / L, where 
k — 1.0 X 10~ 6 in cgs units. The constant C\ specifies the radiative loss through the transition 
region at the footpoint of the half-loop as empirically proportional to the coronal radiative 
loss. AH the constant s used in the model, except ci, in this study are taken the same values 



as m 



Klimchuk et all (120081 ). 



The model assumes that the coronal plasma in a flare loop is heated directly by an 
ad-hoc volumetric heating rate Q(t), or indirectly by non-thermal beams that deposit most 
energy in the lower-atmosphere and then subsequently heat the corona by an upflow (the 
chromospheric evaporation). The M1.0 flare studied here exhibits a small amount of im- 
pulsive X-ray emission above 25 keV energy band for only one minute, so we neglect the 
non-thermal energy input to the two loops in this study. The possible effect by non-thermal 
heating will be discussed in Section HJ The cooling terms include thermal conduction and 
radiation. The model also includes an enthalpy flux, or mass flow. The upward enthalpy 
flux (chromosphere evaporation) plays the role of heating the corona, whereas a downward 
enthalpy flux (coronal condensation) cools down the corona as an energy loss term. The 
enthalp y flux transfers energy and mass between the corona and its base at the transition 
region (IBradshaw fc Cargil]||2010al ). 



Among these terms, the heating rate Q(t) and the loop half-length L are required inputs, 
and the other terms are calculated from the plasma evolution. For L, we assume a semi- 
circular loop, and its half-length is given by L = ttD/A, where D is the distance between the 
tw o footpoints conn ected by the loop. The heating rate Q(t) is inferred using the method 
by iQiu et al.l ( 120121 ). which assumes that the heating flux H(t) = Q(t)L is proportional to 
the UV emission in the rise phase at the footpoint of the flare loop: 



H(t) = A/ exp 



2t 2 



erg cm 2 s , (0 < t < 00) 



(3) 



In the equation, I is the peak count rate (in units of DN s x ) measured by AIA, t is the 
peak time, and r is the characteristic rise time. They are determined by fitting the rise phase 
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of the footpoint UV light curve to a half-Gaussian (see Figure [2]). The heating rate is also 
assumed to be symmetric in time, so it is a full Gaussian even though the fitting uses only 
the rising half of the light curve. The scaling parameter A (in units of erg cm -2 DN _1 ) relates 
the magnitude of the heating flux to the UV peak count rate, and is a free parameter that 
is adjusted to best match the computed plasma coronal radiation to the observed values. 

In this study, we derive H(t) from the four selected footpoints, and therefore model 
four half-loops (LP1, LP2, LP3, and LP4) rooted at the four footpoints. Since these four 
half-loops indeed make two full loops, the parameters L, A, and c\ are the same for the two 
half-loops LP1 and LP4 rooted at FP1 and FP4, and the same for the two half-loops LP2 
and LP3 rooted at FP2 and FP3. The optimal values of A and C\ are chosen from comparison 
with coronal observations. In addition to the impulsive heating rate Q(t), w e also prescribe a 



background heating term Q bk , typically taken as = 10~ 4 erg cm -3 s -1 (IQiu et al.ll2012l ). 
which is a small fraction of Q(t). The EBTEL computes the plasma evolution with the 
initial temperature To and density no- Numerical experiments have shown that selection of 
the initial temperature and density does not modify the plasma evolution once the impulsive 
heating sets in. We use different initial densities in the loops so that the mean plasma velocity 
is zero prior to the impulsive heating, namely, the coronal plasma is in static equilibrium 
before the flare. All these parameters are presented in Table [2j 

Using the inputs described above we compute the coronal response by solving the 
EBTEL equations, (jTJ and 02]), for each half-loop rooted at a UV-brightened footpoint. 
Figure |3] shows the model results. The top four panels show the evolution of the temper- 
ature and density for the four half-loops, and the bottom two panels show the enthalpy 
flux from the corona through its base, with positive flux indicating upflow and negative flux 
indicating downflow. Shown in the figure, the coronal temperature rises with the onset of 
the impulsive heating, and peaks when the heating rate is maximum. Thermal conduction 
front reaching the coronal base then drives the upflow, or chromospheric evaporation, which 
raises the density of the coronal loop. The coronal plasma then cools down as dominated by 
radiative cooling, and in the later cooling phase, the density starts to decrease by coronal 
condensation, or downflow. 

We note that, for each loop system, although we model two half-loops with heating 
functions determined from conjugate footpoints separately, the evolution of the two half-loops 



Table 2: Parameters Used in This Study 



half-loops 


L(Mm) 


A(10 4 erg cm" 2 DN" 1 ) 


Cl 


Q ma x(erg cm" 3 s" 1 ) 


no(10 9 cm 3 ) 


T (MK) 


LP1 & LP4 


28 


1.4 


1.2 


1.6 


0.5 


1.0 


LP2 & LP3 


12 


0.3 


3.5 


1.2 


0.7 


1.0 
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follow along each other very closely. This is due to the fact that the UV light curves at the two 
conjugate footpoints are indeed very similar, with comparable rise time (about 2 minutes) 
and peak count rate. Therefore, the inferred heating rates from conjugate footpoints are 
nearly identical, resulting in very similar evolution of the two half-loops. Such observations 
justify using the EBTEL model for this event. 

On the other hand, the short loop and the long loop exhibit different evolution patterns. 
The short loop is heated for a slightly longer time as derived from the UV rise time, but 
with a smaller heating rate as determined from the best-fit, in combination with a smaller 
length scale that allows a greater thermal conduction. As a result, the peak temperature of 
the short loop is lower than the long loop. The hydrodynamic evolution proceeds faster in 
the short loop, and the turn-over from upflow to downflow takes place earlier in the short 
loop at about five minutes after the onset of the impulsive heating. In the long loop, the 
turn-over takes place at about eight minutes after the heating onset even though the heating 
time is shorter. 



3.2. Comparison with Observations 

Convolving the computed temperature and density with the instrument response func- 
tions, we obtain synthetic coronal radiation fluxes from flare loops to be compared with those 
observed by AIA and EIS. For each full loop, we measure EUV fluxes from a series of 5" x 5" 
boxes along the length of the loop (see Figure [1]). The measured EUV fluxes from these 
different locations exhibit nearly identical time profiles, suggesting that the OD treatment 
by EBTEL is a close approximation in this case. We sum up all the EUV fluxes from these 
boxes as the total emission from the loop, and compare it with the computed EUV flux 
along the full length of the loop. In addition to radiation flux, the model also computes the 
mean enthalpy flux, which is the difference between conduction flux and radiation flux in 
the corona. From the enthalpy flux, we can estimate the plasma flow velocity to compare 
with Doppler velocities measured by EIS at the footpoints of the flare loops. Note that our 
model only has two free parameters, A and c\, which can be determined by best matching 
model results with observed emissions in two AIA channels. In this paper, we choose AIA 
94 A and 335 A for this purpose. Then with this optimal set of parameters, we compare 
model results with observations from other AIA channels and from EIS to verify validity of 
the model. 
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3.2.1. Comparison with AIA Observations 

Figure H] shows the comparison between computed and observed EUV fluxes in AIA 94 
A, 335 A, 211 A, and 171 A. The comparison in 193 A and 131 A bands is not shown due 
to saturation in observed counts in these bands. 

The figure shows that, for both loops, the computed EUV fluxes reasonably agree with 
the observed fluxes in most AIA channels: they all exhibit similar evolution time scales, rise 
and peak around the same time, and their magnitudes are comparable within a factor of 
two. Comparing the evolution of the two loops, it is seen that the short loop evolves more 
quickly, reaching the maximum in higher temperature emission (such as in 94 A) earlier than 
the long loop. The modeled results are basically consistent with the AIA observations. 

There are a couple of discrepancies between model and observations, which are likely 
caused by the limitation of the single temperature approach of the OD model. Early in the 
flare, an impulsive emission peak shows up in nearly all EUV bands. This is produced by 
low-temperat ure (< 1 MK) plasmas a t the footpoint upon impulsive heating at the lower 



atmosphere (IBrosius fe Holmanl |2012| ). In the late phase of the flare loop evolution, the 
observed abundant EUV flux in a few bands (such as in 94 A) is likely produced by cooled 
plasma at lower temperatures of < 1 MK. Since the OD EBTEL model does not repro- 
duce temperature distribution along the loop at a given time, it cannot capture emissions 
at temperatures beyond the averaged value, which are observed by AIA due to its broad 
temperature response. This may also cause deviation of peak time (e.g., in 94 A for the long 
loop). Despite these discrepancies, the overall agreement between model and observations 
suggests that our method reasonably describes mean properties of heating and evolution of 
flare loops. 



3.2.2. Comparison with EIS Observations 

Comparison is also made with EIS observations. Figure [5] shows the synthetic and 
observed emission fluxes in the EIS Fe XV line, which typically forms at 2 MK. It is evident 
that, in both loops, the computed flux tracks the observed flux, with the short loop evolving 
more quickly than the long loop. The absolute magnitude of the observed and computed 
fluxes are comparable by within a factor of two. 

The result for the Fe XII line is not shown. The computed mean temperature is higher 
than the formation temperature of this line, hence it does not produce significant emission 
in this line. In reality, there is a temperature distribution (namely, the differential emission 
measure) in the flare loop, the low temperature portion of which would produce emission in 



- 10 - 



this line. Such cannot be described by the 0D EBTEL model. 

The spectroscopic capability of EIS allows us to further compare the modeled mean 
flow velocity with the Doppler shift measurements to gain insight in the hydrodynamic 
behavior of coronal plasma. Figure [6] shows the Doppler velocities in the Fe XII and Fe XV 
lines measured at the loop footpoints, with the vertical bar marking the range of velocities 
measured along the full loop. The measurements show that, in the impulsive phase, the hot 
Fe XV line is blue-shifted (positive velocity) and the cooler Fe XII line is red-shifted (negative 
velocity), which may be explained by the momentum balance between upflo w (evaporation) 



and d ownflow (condensation) upon impulsive heating of the flare atmosphere (ICanfield et al. 



19901 ). During the decay, both lines are red-shifted, suggesting that coronal condensation 



dominates during the cooling phase. 

EBTEL computes the mean enthalpy flux £ oc Pv, which is the difference between the 
conductive flux and radiative flux in the corona. This mean-property approach cannot repro- 
duce the velocity profiles along the loop and at different temperatures during the impulsive 
heating phase, but can be used to estimate the mean flow velocity in the corona. We assume 
that the coronal loop has a steady-flow, n(T)v(T) = UqVq, and is in dynamic equilibrium 
P(T) = 2k,Bn{T)T = P , where n , P , and vq are EBTEL computed mean density, pres- 
sure, and velocity of the coronal plasma, respectively. This yields the estimate of the flow 
velocity v (T) at a given temperature. 

Figure M shows the flow velocities estimated at the formation temperatures of the Fe 
XII and Fe XV lines in comparison with the measurements. It is seen that the model and 
observations qualitatively agree with each other in the higher temperature Fe xv line, both 
showing the upflow during the impulsive phase and downflow during the decay phase. Specif- 
ically, both model and observations show that the turn-over from upflow to downflow occurs 
around the same time. It occurs after 1:40 UT for the long loop and before 1:40 UT for the 
short loop. 

We note that the magnitudes of the modeled and measured velocities are not directly 
comparable especially during the impulsive phase, when the steady-flow and equilibrium 
assumptions are questionable, and the plasma properties along the loop may largely deviate 
from the mean values computed by EBTEL. Furthermore, the EIS scans missed most part 
of the impulsive phase, when upflow of high temperature plasmas may be predominant. 
Therefore, during the impulsive phase, the comparison is rather poor between the EBTEL 
modeled flow and the observed flow in the two relatively low-temperature lines. 



- 11 - 



4. Summary and Discussions 



We analyze and model two distinctive loop systems that are treated as two big loops in 
an M1.0 flare. Using the empirical heating functions inferred from observed UV light curves, 
we compute the evolution of plasmas in these loops via EBTEL. With the model result we 
reproduce the EUV synthetic fluxes and compare them with observations by AIA and EIS 
for each of the two loops. We also compare the computed mean enthalpy flow velocity with 
the measured Doppler velocity from EIS. The model results agree with observations in many 
aspects. The two different loops show different coronal plasma evolution patterns, with the 
short loop evolving more quickly than the long loop. These patterns are revealed in obser- 
vations and as well reproduced in the model. Our results suggest that the observationally 
constrained heating rates and the EBTEL modeling with only two free parameters provide 
a reasonable description of the mean properties of heating and evolution of flare loops. 

In our model, the loop length, and the timing and duration of the heating rate are derived 
directly from observations. Previous study using the similar method uses exactly the same 
set of parameters to scale the heating rate and energy loss rate for a few thousand loops. Our 
experiment of modeling two spatially resolved loop systems then allows us to investigate the 
variations of these model parameters A and c\ in different loops. The parameter A relates 
the UV emission counts directly from data to the amount of direct heating in the corona 
(QL), and C\ scales the total energy loss through the transition region, or the base of the 
flare loop. Our results show that, to best match the mean properties of loops, different 
parameters have to be used. The best-match A for the long loop is four times that of the 
short loop (see Table [2]). On the other hand, C\ for the long loop is one third that of 
the short loop . For the long l oop, the value of C\ is comparable with the average value 
determined in lOiu et al.l (120121 ) who studied a few thousand long loops; and for the short 
loop, the optimal c\ value is substantially larger. These differences may indicate different 
heating mechanisms and their subsequent low-atmosphere responses in the two loops. A 
plausible explanation for these differences may involve heating by electron beams. Imaging 
observations by RHESSI suggest that during the early impulsive phase, there are hard X- 
ray sources located at the feet of the short loop. In- situ heating by non-thermal electrons 



prod uces significant UV emission at these locations ( lAllred et al.l 120051 ; ICheng et al.l 12010 



20121 . and references therein). Therefore, less amount of UV emission is contributed by direct 
heating in the corona, which heats the lower atmosphere through thermal conduction during 
the impulsive phase. A semi-qu antitative analys is including non-thermal heating term in 
the EBTEL modeling is given by iLiu et al.l ( 120121 ). A larger c\ value used for the short loop 
also indicates a greater amount of radiative loss through the transition region than in long 
loops. This sug gests that Ci might not be an independent constant but may vary with the 
loop dynamics (jCargill et al.ll2012l ). 



-12- 



Flow velocity is an important property in studying the dynamic evolution of flares. 
Spectroscopic velocity measurements can help diagnose the process of energy release and 
plasma evolution. During the impulsive phase, ch romospheric evapo r ation (upflow), driy en 
either by electron beams or thermal conduction (lActon et al.lll982l ; iFisher et al.l ll985aU b; 
Brosius fc Phillipsll2004j ; lMilliganll2008l ). provides indirect heating apart from the direct heat- 
ing (Q) in the corona. In the cooling phase, coronal draining (downflow) dominates. In our 
present study, the spectroscopic observations do not adequately cover the impulsive phase 
of the flare, but our analysis shows that the modeled timing of the turn-over from upflow 
to downflow is consist ent with the measured Dop pler velocity in the decay phase, which 
supports the result of iBradshaw fc Cargill J2010al lbl). Our analysis shows the potential of 
using Doppler shift measurements to constrain the EBTEL model that invokes enthalpy flow 
in the flare evolution. In our future investigation, the high efficiency of the EBTEL model 
will allow us to model a good number of spatially resolved loops, and statistically investigate 
the relationship between energy release and plasma flow. 

Despite its advantages, the EBTEL model does not reproduce the temperature distri- 
bution along the loop, which may explain the insufficient low temperature plasma emis- 
sion from the model, as compared with observations. Furthermore, recent studies suggest 
that the observed evaporation velocity during the impulsive ph ase is temperature dependent 
( IMilligan fc Dennis! 120091 ; IChen &: Ding||2010l ; iLi &: Dingll201ll ). which requires more sophis- 
ticated models, such as ID hydrodynamic model, to compute the plasma evolution especially 
during the impulsive heating phase. The temporal properties of heating rates determined in 
this study can be used as the reference for ID models; on the other hand, a ID model in- 
cludes more assumptions or free parameters such as those describing distribution of heating 
along the loop and therefore invokes addititional observational constraints. 
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Fig. 1. — Top panel: normalized light curves of the 2011 February 16 M1.0 flare in GOES 1-8 
A, AIA 1600 A, and AIA 193 A. The vertical dashed lines denote the time ranges of EIS five 
scans over the flaring region. Middle panels: AIA 1600 A and 193 A images of the flare. The 
red and blue contours in the left panel refer to the positive and negative magnetic polarities, 
respectively. The right panel zooms in the box region of the left panel. The four small black 
boxes mark four footpoints, and the series of small grey boxes are used to measure the total 
emission from the loop. Each box has a cross-sectional area of 5" x 5". Bottom panels: EIS 
Fe XV intensity and Doppler velocity maps for the box region around the GOES 1-8 A peak 
time. 
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Fig. 2. — Observed UV 1600 A light curves (solid) of the four footpoints. The impulsive 
rise phase of each light curve is fitted with a half-Gaussian function. Then the symmetric 
full-Gaussian (dotted) is used to construct the heating function (see text), which corresponds 
to the right coordinate. 
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Fig. 3. — The top four panels show computed evolution of the plasma temperature (solid 
lines) and density (dashed lines) in half-loops rooted at the four footpoints. The bottom two 
panels plot the mean enthalpy flow in these half-loops. 
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Fig. 4. — Comparison between synthetic (dashed lines) and observed (solid lines) EUV fluxes 
of the two full loops in 335 A, 94 A, 211 A, and 171 A, respectively. 
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Fig. 5. — Comparison between synthetic (lines) and observed (triangles) fluxes of the two 
full loops in the Fe XV line. The horizontal bars at triangles represent the time ranges of 
each EIS scan. 
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Fig. 6. — Comparison between computed equivalent velocities from the enthalpy flows (see 
text; solid and dashed lines) and measured Doppler velocities (triangle and asterisk symbols) 
at the four footpoints in the Fe XV and Fe XII spectral lines. The vertical lines at symbols 
represent the observed velocity ranges along the full loop. For clarity, we plot the velocity 
range in the Fe XV line only at FPs 1 and 2, while in the Fe XII line only at FPs 3 and 4. 



